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Abstract 

The effect of randomness or jitter of the acoustic waveform on the 
spectrum of rotor impulsive noise is studied, because of its importance for 
data interpretation. An acoustic waveform train is modelled representing 
rotor impulsive noise. The amplitude, shape, and period between occurances of 
individual pulses are allowed to be randomized assuming normal probability 
distributions. Results, in terms of the standard deviations of the variable 
quantities, are given for the autospectrum as well as special processed 
spectra designed to separate harmonic and broadband rotor noise components. 
Consideration is given to the effect of accuracy in triggering or keying to a 
rotor one per revolution signal. An example is given showing the resultant 
spectral smearing at the high frequencies due to the pulse signal period 
variability. 
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frequency bandwidth, sec - * 

one-sided autospectral function, Eq 2, sec»Pa^ 
see Eqn. 14 

pulse identifier in time records, integer 

/=i 

identifier of data sample, each sample of length T, 
integer 

total number of data samples, integer 

number of rotor periods per sample period, integer 

integers 

sample size, generally equals an integer power of 2 
acoustic pressure time history, Pa 
■ p(nAt), data sequence representing p(t), Eq. 1 
Fourier transform of p(t), sec*Pa 
“ P (f m ) , discrete Fourier transform sequence, 

Eq. 1 

digitization time interval of p(t) data, sec 
» NAt, sample record length, sec 
reference time in sample pressure history 
records , sec 

rotor period of rotation, sec 
standard deviation of variable x 
variation of pulse occurrence time, sec 
phase of P m , rad 
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Subscripts and Superscripts 


J 


i i i 
k m’ km' km 


i 

k T 


value of terms for m-th frequency band for the 
particular i-th pulse of the k-th record sample. 
When term is averaged over k Q samples, "a 11 
replaces "k" . When "i" or "k" is removed, term is 
independent of these • 

variability for i-th pulse of the k-th record 
sample 


Introduction 

The study of helicopter impulsive noise involves analysis of acoustic 
pressure time histories and associated spectra. For a discrete noise source 
such as blade-vortex interaction (BVI), the acoustic pressure pulses in the 
time histories can be randomized somewhat in amplitude and phase with respect 
to the blade passages. This occurs even in well controlled wind tunnel 
environments, where microphones are fixed with respect to model rotor hub 
positions, because of unsteadiness or jitter in blade motion and blade tip 
vortex trajectory. This jitter causes spectral smearing and a broadband 
appearance of the higher harmonics. This can make it difficult to identify in 
an autospectrum the discrete and truly broadband (turbulence caused) noise 
source contributions [1]. In addition, jitter presents a major problem for a 
• common analysis procedure employed to "clean-up" acoustic pressure time 

histories [2]. The procedure involves averaging time histories, in a manner 
which is keyed to the shaft position (such as through a one-per-revolution 
signal), to remove broadband noise resulting in smooth acoustic pressure 
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signatures. Conceptionally spectra corresponding to these averaged signatures 
would represent the "discrete" contribution to the total autospectra. However 
this averaging procedure can render signatures of noticeably distorted shape 
and decreased amplitude unless jitter is small* This concern has prompted 
caution in the selection of test conditions to analyze in this manner. Wisely 
still more caution has been exercised in presentation of special "discrete" 
spectra, being without the benefit of a quantitative understanding of the 
effect of jitter on the spectra. 

It is the purpose of this paper to provide a diagnostic tool by 
quantifying the effect of impulse signal jitter on the narrowband 
auto-spectra. The effect of jitter on special spectral processing, intended 
for separating the discrete from broadband components, is also considered. 

Digital Processing and Terminology 
Digital procedures [3,4] for data processing require that acoustic 
pressure time histories p(t) be sampled at points At apart. The record length 
becomes T * NAt where N is the sample size (typically 1024). The continuous 
p(t) is replaced by the data sequence {p n } “ {p(nAt)} for 
n ■ 1, 2, 3, . . . N. Also replacing a continuous transform P(f) is a 
discrete Fourier transform sequence {P m } « {P(mAf)} ■ {P(f m )} f° r *» “ 1, 

2, 3, ... N. Here Af * 1/T which is the minimum bandwidth resolution. The 
sample size N induces a Nyquist cutoff or folding frequency given by f c * 
l/2At. Prior to digital conversion p(t) is normally low-pass filtered below 
f c to prevent serious aliasing errors. With 211(mAf)nAt - 2 limn /N 
N 

P - P(f ) = At l p exp(-j 2Hmn/T) , m * 1, 2, 3, ... N 
mm . n 

n«l 

“ d N 

p - p(nfit) - if Y P exp( j 2Hm/T) , n - 1, 2, 3, ... N 
n , m 

m-1 


(i) 


which are Fourier transform pairs. It is seen that the data is treated as if 
they were periodic data of period T. From the above P m , the one-sided 
autospectral density function is 

G = G(mAf) - | E [P*Pl 
m T mm 

k k 


_2_ y 
= k T l 
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Here E[ ] is the expected value of the argument and k represents the 
particular data sample each of length T. The total number of averages is 
assumed to be large for the estimate to minimize statistical error. 

Impulsive Noise Transform Decomposition 
Figure 1 illustrates two sample acoustic pressure time histories with four 
pulses in each waveform. This can be representative of a four bladed rotor 
producing BVI, thickness, and/or loading noise. The samples shown are keyed 
to a rotor one-per-revolution signal so each t corresponds to the rotor at a 
particular azimuthal position. Each pulse, identified by the "i" notation, 
can be arbitrary in shape although similar pulses are shown. Presentations 
are also possible for a "b" bladed rotor where the sample period is T = £T r , 
T r being the rotor period of rotation. The number of pulses seen then would 
be b£ rather than 4 for each sample. As shown in Figure 1, variation in 
amplitude and time spacing can occur between the different samples. What is 
desired is to find the spectral consequences of such signal "jitter" between 
samples. 

The first step in the analysis is illustrated in Figure 2 where the 

^p(t) sample of Figure 1 is decomposed into 4 (or b £) separate pulse 
histories, i.e. 

^ < 

. p(nAt) » l , p(nAt) 

K i-1 K 
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The history ^p(nAt) is taken as a pulse signature located at nominal 

(or average) time T with a possible variation of *t. For the moment T 

O it o 

is assumed constant for perfect keying to a one-per-rev signal. The history 

2 2 
. p(nAt) is that pulse at an average time of T + T/4 with variation . x. 

«t O It 

3 4 

The same follows with ^p(nAt) and ^p(nAt) with average spacing T/4 and 
variations ^x. The corresponding transforms are from Equation (1) 


kP “ ’ i-i kP " " t-i 


£ C m - J 2M » (T o + + 


(3) 


which follows from the fact that p(t - t Q ) gives a transform 

P(f )exp[-j2llft 0 ] . Equation (3) explicitly accounts for transform phase 

shifting corresponding to pressure data time shifting. The terms 

fc exp l -j* A ] are Fourier coefficients specifying the shape and 
it id it in 


amplitude of the k-th sample of the i-th pulse. The positive magnitude 


f C ■ |fp I and the phase denoted by f<fc . For the general case of a 

it ID I K m [ it ID 

presentation of "t M periods for a "b" bladed rotor, one replaces the number 4 
in Equation (3) by bi. in the two places appearing. Note then that f m = m/T 
m m/lT r ; so if the number of revolutions analyzed l is say 2 then f m 
equals one-half that for £*1 thereby reducing the bandwidth by a factor of 
two. 


Spectrum for No Jitter Case 

For the limiting case where there is no jitter in amplitude or time for 
the signals, fx * 0 and fc expl-jf* 1 “ *C exp[-j*$ ] for all samples k. It 

lv Iv HI Ev ID ID ID 

is further assumed for this example that the pulses are identical in shape and 

amplitude, i.e. i Cexp[-j i 4] - C exp[-j<(, ] for all pulses i. 
m m m T m 


i 
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From Equation (3) 

bl 

l 

i«l 


P - I C exp[-j<fr - j2wf (T + )] 

m m Kl JT m J m o bl J 


b£C exp[-j$ - j2irf T ] for m * b£, 2b£, 3b£, 
m m mo 


for m * b£, 2b£, 3b£, . . . 


It is seen that a nonzero contribution is obtained only for frequencies that 
are multiples of f^> where f^ * b£/T * b/T^ which is the blade passage 


frequency. The auto-spectrum corresponding to P m is, from Equation (2), 

for m * b£, 2b£, 3b£, ... 


G . - v 2 

■ 0 for m * bJt, 2b£, 3b£, ... 

Spectrum for a Normal Proability Distribution of Jitter 
Now considering the more general case where jitter is present, from 
Equation (2) and (3) 


(4) 


2 1 r * 

— — ) P P 

T k 4 km k n 
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This equation for 6 is now evaluated by assuming pulse amplitude, shape, 

m 

* 

and time of occurrence are all random variables. The term [, P ,P ] is 

k m k m 

a function of *C , , and fx. For the moment call these 

fc m ic m K 

55 j, ^2 » and respectively, each representing b£ random variables (for 

example Xj - 1^, k C m ,.., k c m D* 

For large ko 


Q ^ CO OO oo ^ 

r— 7 [.P ,P]“ / / J[PP] p(x,, x 0 , x_)dx. dx 0 dx_ 
k , L , k m k m 1 1 J m m ~1 *2 *3 *1 ~2 

o k a l — 00 — oo — * 


( 6 ) 


8 


The p(x . x , x„) term is the joint probability density. It is now assumed 
* a/ 1 " a^2 ~3 

that x , x , and x_ variations are uncorrelated and additionally that their 

/w 1 ’ ^ 2 . 7 

variations all have normal or Gaussian distributions, i.e. 
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= n 

i=l 
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■ n 
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P( k T> 


where for each component of x^, x^, and x^ 

p(x) = — exp[-(x - x ) 2 /2 a 2 ] . 

a /H ox 

x 


(7) 


( 8 ) 


Here cr x is the standard deviation of x and Xq is the mean value of x, i.e. 


x = E[x] - / x p(x) dx 

o 


(9) 


and 

a 2 = E [(x - x ) 2 ] “ E [x 2 ] - x 2 . (10) 

X o o 

It is acknowledged that the properties assumed in Equations (7) and (8) are 
restrictive in allowable pulse behavior. It does, however, permit a solution 
believed to be realistic in practice. 

.It is seen from Equations (6) and (7) that G m of Equation (5) can be 
evaluated by the individual determinations of the expected values of the 
separate terms. First employing Equations (9) and (10) 


i i» 12 2 

E [ C C ] - C * o. , if i - i* 
m m am i c ’ 
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(ID 


i i 1 

C C , if i * i' 

a m a m 


i i 1 i 

where X C ■ E [ C ] ■ ■; — T . C which is the average amplitude of the 
a in m lc . . tc in 

o k=l 


coefficient and the o-term is the standard deviation of C . 

m 
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Next for the phase 


E [exp[-j( i i - V )]] - 1 if i - 1' 
m m 


exp[-j ( - ^ m )] exp [- y (o.? + a ± 2 )] 


if i * i' 


m 


m 


where $ is the average phase of the coefficient and the a-term its 
a m 


( 12 ) 


standard deviation. Finally for the time jitter term 
E [exp(-j2llf m ( - i x))] - 1 if i - i’ 

- exp [- y(2JIf m ) 2 ( o 2 + o 1 2 )] 
if i * i'. T T 

Now inserting the results of Equations (11) - (13) into Equation (5) and 

rearranging one obtains 

G - (G ) + (G ) 

m m r m d 


(13) 


(14) 


where 


, bt . , 

(G) r I V [1 + (o VC t) - exp[-o " - (2Hf a. )*]] 

m r i , , a m a m i, ml 


2 ,i„2 


and 


i-1 


bt bt 
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r m 


(G ). - ~ l l i C i C exp[-j( - % ) - j2H ( ] 

m d T i*»l a m am v J a r m a T m J bt 

• exp [-y (o v 2 + a ± 2 ) -y (2nf ffl ) 2 (a 2 , + a ± 2 )] . 


m 


m 


This is the resultant autospectrum where the statistics of jitter are defined 

in terms of the standard deviation of the variable quantities - pulse time, 

amplitude, and shape. The term (G,,,)^, defines a so-called "discrete" 

portion whereas (G ) defines a "random" portion of G . This is dealt with 

m r m 


below 
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Special Processing 

In the proceeding analyses T 0 of Equation (3) is assumed constant for 

clarity of presentation. However, even if data samples were not perfectly (or 

at all) keyed to a rotor one-per-rev signal through say triggering error and 

thus T would be variable between samples k, i.e. T = T in Equation 
o o k. o 

(3), it is readily shown that Equation (14) for is unchanged. 

Such keying, however, would have impact on certain alternate spectral 
definitions one might employ to study particular noise mechanisms. Of 
interest is the spectra corresponding to averaged time histories. As 
indicated in the introduction, such spectra may be produced with the intent to 
remove broadband noise content. This "discrete" or deterministic 

f 

spectrum, call this (G ) , , is found to be 

m a 


(G ). 
m a 



2 


“ (G ) exp [-(2llf a _) 2 ] 
m a mi 

o 

where a T is the standard deviation of T (assuming normal and 
o o 

independent variability) and (G ), is that from Equation (14). In 

m a 

parallel one can define a "random" spectrum (G^)^ equal 


(15) 


(G )’ - | E P - E[P ] 
m r T m m 


G - (G ) 
m m a 


(16) 


where G is that from Equation (14). So only when variations in T are 
m o 


minimized do (G )' * (G ), and (G ) f * (G ) • It is seen that 
md md mr mr 

Equations (15) and (16) relate the degree of success one can have in 

processing to separate ,, discrete n from "random" noise, to the extent such 
success depends on keying accuracy. Equation (14) can be employed to help 
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determine to what degree one can utilize the results, even with accurate 
keying, for any physical interpretation of the separation spectral data. 

Example 

The spectrum for a signal with jitter is now compared to that spectrum 
obtained for an ideal case where there is no jitter. For the no jitter case, 
all o's are zero and if the pulses are identical for different i's then 
Equation (14) renders Equation (4). For the comparison it is assumed for the 
jitter case that all pulses are of the same average shape, i.e. 


*C exp[-j*<j> ] * C exp[-j <j> ). Additionally, it is assumed any 
dm dm dm dm 

amplitude vdriations are uniform for each pulse so no shape charges occur, 


only overall amplitudes, i.e. a. ,i is constant (say a c /c) and o - 0. 

1 G / c - 


m 


m 


Also the time variability for the different pulses i are the same, so 


0^*0^. So for the jitter case from Equation (14) 

T 1 


G m 53 : >, + < G > < 17 > 

m m a m r 

(G ) = bl C 2 [l + (— ) Z - exp [ — ( 2IIf a ) Z ]] 
mr T am c r rat 

and 


(G ) . - (b£ C ) 2 exp [~(2IIf a ) 2 ] for m - b£, 2 bi, 3 b£, . . . 

m a l a m m x 

=0 for m * b£, 2 b£, 3 b£, . . . 

These spectra and that of Equation (4) are presented in Figure 3 for the first 

m = 34 frequency bands. The levels are normalized in each band by 

2 

2(b£ C ) /T which correspond, at harmonics of f, , to the levels of G 
am Dm 

for the ideal case of no jitter (Equation (14)). The data period T is taken 

as equal to the rotor period T r (£“1), the number of rotor blades are four 

(b=4), and the time variability is large at <j t /T » 0.01. For clarity. 



12 


the amplitude variability a c is assumed zero so only the effect of o T is 
shown. In general the a contribution to the comparison would be 
2 

10 log [(a c /c) /bfc]. For the present b£.“4 case, a value of oj c=0.5 
would render a -12dB contribution for all m bands. 

Figure 3 shows the effect of pulse time jitter to be significant for the 
spectrum G m . The levels of the individual blade passage harmonics are 
substantially reduced and the spectra becomes broadband in appearance with 
levels approaching -6 dB (-10 log b£) below the no jitter case for large m. 
Note in this normalized presentation that the integrated levels for the 
spectra with and without jitter are equal, promoting the view that jitter 
merely shifts the spectral energy content from the harmonics of fb to the 
side bands. However, this does not strictly hold, for as indicated in 
Equation (17) for the unnormalized G m , the level of each band m is governed 
by the coefficient amplitude a C m for that band, not that of the nearest 
band corresponding to a harmonic of fb« 

Conclusions 

The present results explain quantitatively the spectral smearing which 
occurs when simple jitter is present in the acoustic signals. The results can 
be applied diagnostically in the identification and separation of source 
contributions to the rotor noise spectra. Applicability of the present 
equations is direct when discrete or impulsive noise dominates the noise 
field. In the presence of additional discrete noise with differing pulse 
shape and jitter statistics, as well as the presence of substantial broadband 
noise, a precise analysis requires additional steps. Here appropriate terms 
can be readly added to the modelled pressure transform of Equation (3) 
prior to determining the autospectrum or special spectra processed for 
diagnostic purposes. 
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